February 2, 2010 16:20 WSPC/INSTRUCTION FILE paper 



International Journal of Modern Physics E 
© World Scientific Publishing Company 



MEDIUM EFFECTS ON CHARGED PION RATIO IN HEAVY ION 

COLLISIONS 



CHE MING KO 

Cyclotron Institute and Department of Physics and Astronomy, Texas A&M University, College 

Station, Texas 77843-3366 
ko&comp. tamu. edu 

YONGSEOK OH 

Korea Institute of Science and Technology Information, Daejeon 305-806, Korea 

yoh@kisti.re.kr 

JUN XU 

Cyclotron Institute, Texas A&M University, College Station, Texas 77843-3366 
xujun@comp.tamu.edu 

Received (received date) 
Revised (revised date) 

We have recently studied in the delta-resonance-nuclcon-holc model the dependence of 
the pion spectral function in hot dense asymmetric nuclear matter on the charge of the 
pion due to the pion p-wave interaction in nuclear medium. In a thermal model, this 
isospin-dependent effect enhances the ratio of negatively charged to positively charged 
pions in neutron-rich nuclear matter, and the effect is comparable to that due to the 
uncertainties in the theoretically predicted stiffness of nuclear symmetry energy at high 
densities. This effect is, however, reversed if we also take into account the s-wave inter- 
action of the pion in nuclear medium as given by chiral perturbation theory, resulting 
instead in a slightly reduced ratio of negatively charged to positively charged pions. Rel- 
evance of our results to the determination of the nuclear symmetry energy from the ratio 
of negatively to positively charged pions produced in heavy ion collisions is discussed. 



1. Introduction 

The nuclear symmetry energy is the energy needed per nucleon to convert all pro- 
tons in a symmetric nuclear matter to neutrons. Knowledge on the density depen- 
dence of nuclear symmetry energy is important for understanding the dynamics of 
heavy ion collisions induced by radioactive beams, the structure of exotic nuclei 
with la rge neut ron or proton excess, and many important issues in nuclear astro- 
physics l l | 2P | 4 | At normal nuclear matter density, the nuclear symmetry energy 
has long been known to have a value of about 30 MeV from fitting the binding 
energies of atomic nuclei with the liquid-drop mass formula. Somewhat stringent 
constraints on the nuclear symmetry energy below the normal nuclear density have 
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also been obtained during past few years from studies of the isospin diffusion 
and isoscaling^in heavy-ion reactions, the size of neutron skins in heavy nuclei j^j , 
and the isotope dependence of giant monopole resonances in even- A Sn isotopes L^M. 
For nuclear symmetry energy at high densities, transport model studies have shown 
that the ratio of negatively to positively charged pions p roduce d in heavy ion colli- 
sions with neutron-rich nuclei is sensitive to its stiffness I 12 | 13 | Comparison of this 
ratio from an isospin-dependent Boltzmann-Uchling-Ulenbeck (IBUU) transport 
model base d on the non-relativistic momentum-dependent (MDI) nuclear effective 
interactions with measured data from heavy ion collisions by the FOPI Col- 
laboration a t Qgj seems to indicate that the nuclear symmetry energy at high 
density might be very soft. Although this study does not include the relativistic 
effects, which may affect the charged pion ratio as shown in Ref it provides 
an important step in the determination of the nuclear symmetry energy at high 
densities. 

The transport model used in Ref. ^1 neglects, however, medium effects on pi- 
ons, although it includes those on nucleons and produced A resonances through 
their isospin-dependent mean- field potentials and scattering cross sections. It is 
well-known that pions interact strongly in nuclear medium as a result of their 
p-wave couplings to the nucleon-particle-nucleon-holc and dclta-particlc-nuclcon- 
hole (A-hole) excitations, leading to the softening of their disp ersion relations or 
an increased strength of their spectral functions at low energies | 17 | 18 | 19 | 20 | 21 | 22 ] 
Including pion medium effects in the transport model has previously been shown 
to enhance the production of low energy pions in h igh energy heavy ion collisions, 
although it does not affect the total pion yield ^31. Since pions of different charges 
are modified diffe rent ly in asymmetric nuclear matter that has unequal proton and 
neutron fractions including such isospin-dependent medium effects is expected 
to affect the ratio of negatively to positively charged pions produced in heavy ion 
collisions. 



2. Pion p-wave interactions in nuclear medium 

Considering only the dominant A-hole excitations as in Ref. ^51, as the contribution 
from the nucleon particle-hole excitations is known to be small, the sclf-cnergy of a 
pion of isospin state mt, energy w, and momentum A: in a hot nuclear medium due 
to its p-wave interaction is given by 

K^U^X^k) £ \(lm T \lm t lm T )\ 2 

^ 77 ' m T ,m T 

d 3 p 1 ( 1 1 \ 



(1) 
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' m A + U A T + (k± p) 2 /2m A - iT A T /2 - m N - - p 2 /2m N . In 
the above, m T ~ 138 MeV, m N ~ 939 MeV, and m A ~ 1232 MeV are the masses 
of pion, nucleon, and A resonance, respectively; / A ~ 3.5 is the ttNA coupling 
constant and F 7r (k) = [1 + 0.6(fc 2 /m^)] -1 / 2 ^ is the 7riVA form factor determined 
by fitting the decay width Ta — 118 MeV of A resonance in free space. The 
summation in Eq. (fTJ) is over the nucleon isospin state m T , and the A resonance 
isospin state m T ; and the factor (| m T \l mt \ m T ) is the Clebsch-Gordan coefficient 
from the isospin coupling of pion with nucleon and A resonance. The momentum 
integration is over that of nucleons in the nuclear matter given by a Fermi-Dirac 
distribution with p B and pq being, respectively, the baryon and charge chemical 
potentials determined by charge and baryon number conservations; p^ T and U^ T 
are, respectively, the density and mean-field potential of nucleons of isospin state 
m T in asymmetric nuclear matter; and r A T and U A T are, respectively, the width 
and mean-field potential of A resonance of isospin state m T . 

For the nucleon mean-field potential U^ T , we have used the one obtained from 
the momentum-independent (MID) interaction ^ i.e., U^ T {p B ,8\\ kc ) = ot(p B /p ) + 
KPb/Po) 1 + U^{p B ,5 m «), with U^(p Bl 5 hkc ) = -4{F(x)(p B /p ) + [18.6 - 
F{x)]{p B / Po ) G ^}m T 6 likc + [18.6 - F(x)][G{x) - l]{p B / 'p ) G ^S ly J being the nu- 
cleon symmetry potential. The parameters a = —293.4 MeV, (3 = 240.1 MeV, and 
7 = 1.216 are chosen to give a compressibility of 212 MeV and a binding energy 
per nucleon of —16 MeV for symmetric nuclear matter at the saturation or normal 
nuclear density p = 0.16 fm~ . The nucleon symmetry potential U™y~ (p B , Sy ikc ) 
depends on the baryon density p B = p n + p p + p A _ + p Aa + p A+ + p A++ and the 
isospin asymmetry (5 likc = (p n - p p + p A _ - p A++ + p A0 /3 - p A +/3)/p B of the 
asymmetric hadronic matter, which is a generalization of the isospin asymmetry 
^ = (Pn ~ Pp)/i.Pn + Pp) usually defined for asymmetric nuclear matter without 
A resonances ^21. The nucleon mean-field potential also depends on the stiffness 
of nuclear symmetry energy through the parameter x via the functions F(x) and 
G(x). We consider the three cases of x = 0, x = 0.5, and x = 1 with correspond- 
ing values F(x = 0) = 129.98 and G(x = 0) = 1.059, F(x = 0.5) = 85.54 and 
G(x = 0.5) = 1.212, and F(x = 1) = 107.23 and G(x = 1) = 1.246. The resulting 
nuclear symmetry energy becomes increasingly softer as the value of x increases, 
with x = 1 giving a nuclear symmetry energy that becomes negative at about 
3 times the normal nuclear matter density. These symmetry energies reflect the 
uncertainties in the theoretical predictions on the stiffness of nuclear symmetry en- 
ergy at high densities. For the mean-field potentials of A resonances, their isoscalar 
potentials are assumed to be the same as those of nucleons, and their symmetry 
potentials are taken to be the average of those for neutrons a nd p rotons with weight- 
ing factors depending on the charge state of A resonance^, i.e., U A y = U£ sy , 
U A+ = -U p 4- -U n U A " = -U p + -U n and U A ~ = U n 

^ asy 3 ^ asy ' ^ as y ' as y 3 as y ^3 asy ? ^ asy ^ asy ■ 

Including the short-range A- hol e repulsive interac tion via the Migdal parameter 
g', which has values 1/3 < g' < 0.6 1 17 | 18 | 19 | 20 | 21 | 22 ] moc iifi es the pion self-energy 
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to n mt = n™7(l- 3 'II™7fc 2 ). The pion spectral function 5™' (u, k) is then related 
to the imaginary part of its in-medium propagator D mt (u>, k) = l/[w 2 — k 2 — m 2 — 
n mt (w,fc)] viaS; nt (u;,fc) = -(l/7r)ImL> mt (a;,fc). 

The modiheation of the pion properties in nuclear medium affects the decay 
width and mass distribution of A resonance. For a A resonance of isospin st ate 
and mass M and at rest in nuclear matter, its decay width is then given by^S 

(2) 



. nit 



xlm 



\ e (u)-m t u Q )/T _ ^ j y g 

k 2 D mt (cj,k) /2 n mt (cj,fc) 



(m N +k 2 /2m N +U™ T -n B -2m T fi Q )/T , ^ 



3 (l-g'U^(^k)/k 2 ) 



9 



In the above, the first term in the last line is due to the decay of the A resonance 
to pion but corrected by the contact interaction at the ttNA vertex, while the 
second term contains the contribution from its decay to the A-hole state without 
coupling to pion. The first two factors in the momentum integral take into account, 
respectively, the Bose enhancement for the pion and the Pauli blocking of the nu- 
cleon. To include possible chemical non-equilibrium effect, a fugacity parameter z n 
is introduced for pions. The pion energy uj is determined from energy conservation, 
i.e., M + U2 T = w + m N + k 2 /2m N + U^ T . The resulting mass distribution of A 
resonances is then given by P A {M) = A[r^ T (Af)/2]/[(A7 - m A ) 2 + r^ T2 (M)/4], 
where A is a normalization constant to ensure the integration of Pa(M) over M is 
one. 

We have solved Eqs. (Q]) and ([2]) sclf-consistcntly to obtain the pion spectral 
functions and the mass distributions of A resonances in asymmetric nuclear matter. 
The results obtained with the Migdal parameter g 1 = 1/3 are illustrated in Fig. Q] 
and Fig. [2] for an asymmetric nuclear matter of isospin asymmetry <5nkc — 0.133, 
twice the normal nuclear matter density p B = 2p Q , temperature T ~ 43.6 McV, 
and chemical potentials /i B ~ 941.89 MeV and /Iq ~ —18.26 MeV, corresponding 
to those to be used in our thermal model and also similar to those reached in 
the transport model with the nuclear sy mme try energy x = 1 for central Au+Au 
collisions at the beam energy of 0.4 AGevQSl. Shown in Fig.Q] are the pion spectral 
functions as functions of pion energy for different values of pion momentum. It is 
seen that for low pion momenta the spectral function at low energies has a larger 
strength for 7r _ (dotted line) than for 7r° (solid line), which has a strength larger 
than that for 7r + (dashed line). This behavior is reversed for high pion energies. 
Fig. [2] shows the mass distributions of A resonances at rest in asymmetric nuclear 
matter as functions of mass. One sees that they are similar to that in free space 
(solid line) as a result of the cancelation between the pion in-medium effects, which 
enhance the strength at low masses, and the Pauli-blocking of the nuclcon from delta 
decay, which reduces the strength at low masses. This is consistent with the observed 
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Fig. 1. (Color online) Spectral functions of pions in asymmetric nuclear matter of density 2p and 
isospin asymmetry <5i;i; = 0.133 as functions of pion energy for different pion momenta of (a) m^, 
(b) 2m n , (c) 3m n , and (d) 4m n . All are calculated with the Migdal parameter g' = 1/3. 



similar energy dependence of the photo -pro ton and photo-nucleus absorption cross 
sections around the A resonance mass Furthermore, the strength around the 
peak and near the threshold of the A resonance mass distribution slightly decreases 
with increasing charge of the A resonance due to nonzero isospin asymmetry of the 
nuclear medium. 



3. Charged pion ratio in hot dense asymmetric nuclear matter 

To see the above isospin-dependent pion in-medium effects on the n~ /tt + ratio in 
heavy ion collisions, we have used a thermal model w hich assumes that pions are 
in thermal equilibrium with nucleons and A resonances ^1 In terms of the spectral 
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Fig. 2. (Color online) Mass distributions of A resonances at rest in asymmetric nuclear matter of 
density 2p and isospin asymmetry &uk c = 0.133. The solid line corresponds to that in free space. 
The distributions near the threshold and at the peak are enlarged in the insets. 

function Si(u, k), the density of a particle species i is then given by 

P ^ 9 * l (^k^fo ^ z -^-B„!-Q^)/T ± 1 ■ ( 3 ) 

In the above, g iy S,, and Qi are the degeneracy, baryon number, and charge of 
the particle. The fugacity parameter z i is introduced to take into account possible 
chemical non-equilibrium effect. The exponent n i is 2 for pions and 1 for nucleons 
and A resonances. For the spectral functions of A resonances, we neglect their 
momentum dependence and replace the integration over the energy uj by that over 
mass. The uj in the Fermi-Dirac distribution for A resonances is then simply u) — 
AI + k 2 /2M + U^ T . For nucleons, their spectral functions are taken to be delta 
functions if we neglect the imaginary part of their self-energies, i.e., S^ T (u!,k) = 
8(w - m N - k 2 /2m N - U^ T ). 

According to studies based on the transport model USEES ^ c total number of 
pions and A resonances in heavy ion collisions reaches a maximum value when the 
colliding matter achieves the maximum density, and remains essentially constant 
during the expansion of the matter. For Au+Au collisions at the beam energy of 
0.4 AGeV, for which the n~ / tt + ratio has been measured by the FOPI Collaboration 
at GSI L^, the IBUU transport model gives a maximum density that is about 
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twice the normal nuclear matter density and is insensitive to the stiffness of the 
nuclear symmetry energ y, as it is mainly determined by the isoscalar part of the 
nuclear equation of state This density is thus used in the thermal model. The 
temperature in the thermal model is determined by fitting the measured pion to 
nuclcon ratio, w hich is about 0.014 including pions and nuclcons from the decay 
of A resonances without medium effects and with unity fugacity parameters 
for all particles, and the value is T ~ 43.6 MeV. The assumption that pions and 
A resonances are in chemical equilibrium is consistent with the short chemical 
equilibration times estimated from the pion and A resonance production rates. The 
isospin asymmetry of the hadronic matter is then taken to be <5nke — 0.080, 0.106, 
and 0.143, corresponding to net charge densities of 0.920p , 0.894p and 0.85 7,00, 
for the three symmetry energies given by x — 0, 0.5, and 1, respectively, in order to 
reproduce the ir~ /ir + ratios of 2.20, 2.40, and 2.60 predicted by the IBUU transport 
model of Ref. using corresponding symmetry energy parameters without pion 
in-medium effects. Since the medium effects enhance the pion and A resonance 
densities, to maintain the same pion to nucleon ratio as the measured one requires 
the fugacity parameters for pions and A resonances to be less than one. Also, the 
pion in-medium effect s h ave been shown to affect only slightly the pion and the A 
resonance abundance ^3, indicating that both pions and A resonances arc out of 
chemical equilibrium with nucleons when medium effects arc included, as expected 
from the estimated increasing pion and A resonance chemical equilibration times 
as a result of the medium effects. Because of the small number of pions (about 
0.3%) and A resonances (about 1.1%) in the matter, the density, temperature, 
and net charge density of the hadronic matter are expected to remain unchanged 
when the pion in-medium effects are introduced. They lead to, however, a slight 
reduction of the isospin asymmetry to <5nk c — 0.073, 0.097, and 0.133 for the three 
symmetry energies, respectively. We note that with the fugacity of nucleons kept 
at zpf = 1, the fugacity parameters of about z^ = 0.061 and z A = 0.373 are needed 
to maintain same total number of pions and A resonances as in the case without 
pion in-medium effects, and that the required values for the fugacity parameters 
increase only slightly for the other two symmetry energies considered here. 

Results on the ir~ /ir + ratio in Au+Au collisions at the beam energy of 0.4 AGeV 
are shown in Fig. [3] With the value g' = 1/3 for the Migdal parameter, values for 
the 7T~ /ir + ratio are 2.32, 2.60, and 2.94 for the three symmetry energy parameters 
x = 0, 0.5, and 1, respectively, which are larger than corresponding values for the 
case without including the pion in-medium effects as shown by those for g' = oo 
in Fig. [3] These results indicate that the isospin-dependent pion in-medium effects 
on the charged pion ratio are comparable to those due to the uncertainties in 
the theoretically predicted stiffness of the nuclear symmetry energy. The measured 
7r~/7r + ratio of about 3 by the FOPI Collaboration, shown in Fig. [3] by the dash- 
dotted line together with the error bar, which without the pion in-medium effects 
favors a nuclear symmetry energy softer than the one given by x = 1, is now best 
described by a less softer one. 
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Fig. 3. (Color online) The 7r~ /71-+ ratio in Au+Au collisions at the beam energy of 0.4 AGeV 
for different values of nuclear symmetry energy (x = 0, 0.5, and 1) and the Migdal parameter 
g' = 1/3, 0.4, 0.5, and 0.6 in the A-hole model for the pion p-wave interaction. Results for g' = oo 
correspond to the case without the pion in-medium effects. 



Fig. |3] further shows the results obtained with larger values of g' = 0.4, 0.5 
and 0.6 for the Migdal parameter. It is seen that the isospin-dependent pion in- 
medium effects are reduced in these cases compared to the case of g' = 1/3 as the 
repulsive interaction between A-hole states becomes stronger, thus reducing the 
pion in-medium effects. With these larger values of g' , symmetry energies softer 
than that given by x — 1 are then needed to describe the measured n~ /n + ratio. 



4. Pion s-wave interactions in nuclear medium 

The above study does not include the s-wave interactions of pions with nucleons. 
Calculations based on the chiral perturbation theory have shown that the pion s- 
wave interaction modifies the mass of a pion in nuclear medium, and for asymmetric 
nuclear matter this effect depends on the charge of the pion^l Up to the two- loop 
approximation in chiral perturbation theory 1^21 the self energies of n~ , 7r + , and 
7T° in asymmetric nuclear matter of proton density p p and neutron density p n are 



2, 2010 16:20 WSPC/INSTRUCTION FILE paper 



Medium effects on charged pion ratio in heavy ion collisions 9 

given, respectively, by 

n ( Pp , Pn ) = p n [T~ N - T+ N ] - p p [T~ N + T+ N ] + n rcl (p p , p n ) + IL7 or (p p , Pn ) 
H + {p P ,Pn) = TT(p n ,p p ) 

n° ( Pp , p„) = - (Pp + p n )T+ N + n° or ( Pp ,p n ). (4) 

In the above, T ± are the isospin-even and isospin-odd 7riV-scattering T-matrices 
which have the empirical values T~ N w 1.847 fm and w —0.045 fm extracted 
from the energy shift and width of the Is level in pionic hydrogen atom. The term 
IT~i is due to the relativistic correction, whereas the terms 11^,. and iI° or are the 
contributions from the two- loop order in chiral perturbation theory. Numerically, 
it was found in Ref. ^1 that changes of pion masses in asymmetric nuclear matter 
of density p = 0.165 fm -3 and isospin asymmetry S = 0.2 are Am,- = 13.8 MeV, 
Am 7r+ = -1.2 MeV, and Am,. = 6.1 MeV. 
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Fig. 4. (Color online) Similar to Fig. [3] with both pion s-wave and p-wave interactions included. 



Taking into account the isospin-dependent pion self energies due to pion s-wave 
interactions in asymmetric nuclear matter changes the results shown in Figs. [T] and 
__J For the pion spectral function, the one for ir + now has a larger strength at low 
energies than that for 7r~. Similarly, the strength near the threshold of the A res- 
onance mass distribution now increases with increasing charge of the A resonance, 
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although that around the peak still decreases with increasing A resonance charge. 
As a result, the ir~ /ir + ratio in Au+Au collisions at the beam energy of 0.4 AGcV 
is slightly reduced after the inclusion of both pion s-wave and p-wave interactions 
in asymmetric nuclear matter as shown in Fig. [4] 

5. Summary 

The pion spectral function in asymmetric nuclear matter becomes dependent on the 
charge of a pion. For the p-wave interaction of the pion, modeled by its couplings to 
the A-hole excitations in nuclear medium, it leads to an increased strength of the 
7T~ spectral function at low energies relative to that of the tt + spectral function in 
dense asymmetric nuclear matter. In a thermal model, this isospin-dependent effect 
increases the tt~ /ir + ratio from heavy ion collisions, and the effect is comparable 
to that due to the uncertainties in the theoretically predicted stiffness of the nu- 
clear symmetry energy at high densities. However, including also the pion s-wave 
interaction based on results from the chiral perturbation theory reverses the isospin- 
dependent pion in-medium effects, leading instead to a slightly reduced n~ ratio 
in neutron-rich nuclear matter. Taking into consideration of the isospin-dependent 
pion in-medium effects in the transport model thus would have some influence on 
the extraction of the nuclear symmetry energy from measured ir~ /ir + ratio. 
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